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Abstract. We theoretically investigate the localization of an expanding Bose- 
Einstein condensate with repulsive atom-atom interactions in a disordered 
potential. We focus on the regime where the initial inter-atomic interactions 
dominate over the kinetic energy and the disorder. At equilibrium in a trapping 
potential and for the considered small disorder, the condensate shows a Thomas- 
Fermi shape modified by the disorder. When the condensate is released from 
the trap, a strong suppression of the expansion is obtained in contrast to the 
situation in a periodic potential with similar characteristics. This effect crucially 
depends on both the momentum distribution of the expanding BEC and the 
strength of the disorder. For strong disorder as in the experiments reported 
by D. Clement et ai, Phys. Rev. Lett. 95, 170409 (2005) and C. Fort et 
al., Phys. Rev. Lett. 95, 170410 (2005), the suppression of the expansion 
results from the fragmentation of the core of the condensate and from classical 
reflections from large modulations of the disordered potential in the tails of the 
condensate. We identify the corresponding disorder-induced trapping scenario for 
which large atom-atom interactions and strong reflections from single modulations 
of the disordered potential play central roles. For weak disorder, the suppression 
of the expansion signals the onset of Anderson localization, which is due to 
multiple scattering from the modulations of the disordered potential. We compute 
analytically the localized density profile of the condensate and show that the 
localization crucially depends on the correlation function of the disorder. In 
particular, for speckle potentials the long-range correlations induce an effective 
mobility edge in ID finite systems. Numerical calculations performed in the mean- 
field approximation support our analysis for both strong and weak disorder. 
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1. Introduction 



1.1. Disorder and ultracold atomic gases 

Understanding the effect of disorder in physical systems is of fundamental importance 
in various domains, such as mechanics, wave physics, solid-state physics, quantum 
fluid physics, or atomic physics. Although in many situations this effect is weak 
and can be ignored in first approximation, it is not always so. Strikingly enough, 
even arbitrarily weak disorder can dramatically change the properties of physical 
systems and result in a variety of non-intuitive phenomena. Many of them are not 
yet fully understood. Examples in classical systems include Brownian motion [JJ, 
percolation [2], and magnetism in dirty spin systems [31 El [SI 151 HI- In quantum 
systems the effects of disorder can be particularly strong owing to the complicated 
interplay of interference, particle-particle interactions and disorder. The paradigmatic 
example is (strong) Anderson localization of non- interacting particles [SI [10l Qj] . 
Other interesting effects of disorder in quantum systems include weak localization 
and coherent back-scattering [12) , disorder-driven quantum phase transitions and the 
corresponding Bose glass [131 HH US] and spin glass [TH [J7] phases. 

Anderson localization (AL) signals out in two equivalent ways, either as the 
suppression of the transport of matterwaves in disordered media, or as an exponential 
decay at large distances of the envelope of the eigenstates of free-particles in a 
disordered potential [IT] . Both properties strongly contrast with the case of periodic 
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potentials, in which transport is free and all eigenstates extend over the full system as 
demonstrated by the Bloch theorem [T5] . Anderson localization is due to a destructive 
interference of particles (waves) which multiply scatter from the modulations of a 
disordered potential. It is thus expected to occur when interference plays a central 
role in the multiple scattering process [11] . In three dimensions, it requires the particle 
wavelength to be larger than the scattering mean free path I as pointed out by Ioffe 
and Regel [19]. One then finds a mobility edge at momentum k — l/l, below which 
AL can appear. In one and two dimensions all single-particle quantum states are 
predicted to be localized [50] [5T] [55] , although for certain types of disorder one has an 
effective mobility edge in the Born approximation [23, 241 125j. 

Ultracold atomic gases are now widely considered to revisit standard problems of 
condensed matter physics under unique control possibilities. Dilute atomic Bose- 
Einstein condensates (BEC) [26j [27] EEl [29] and degenerate Fermi gases (DFG) 
[5UI I5TI 1551 1551 154] are produced routinely taking advantage of the recent progress in 
cooling and trapping of neutral atoms [35] [33 [37] ■ In addition, controlled potentials 
with no defects, for instance periodic potentials (optical lattices), can be designed 
in a large variety of geometries [58] . In periodic optical lattices, transport has been 
widely investigated, showing lattice-induced reduction of mobility (39] [40] [4Tj and 
interaction-induced self-trapping (42] [43] . Controlled disordered potentials can also be 
produced optically as demonstrated in several recent experiments [HI 2S] [35] H7J SH] , 
for instance using speckle patterns [49] [50]. Other techniques can be employed to 
produce controlled disorder such as the use of magnetic traps designed on atomic 
chips with rough wires [SH E2] [53] [54 , 55J, the use of localized impurity atoms [56"[l57]. 
or the use of radio- frequency fields [58]. However, the use of speckle potentials has 
unprecedented advantages from both practical and fundamental points of view. First, 
they are created using simple optical devices and their statistical properties are very 
well known [59] [60]. Second, they have finite-range correlations which offers richer 
situations than theoretical (^-correlated potentials (i.e. uncorrelated disorder) and 
the correlation functions can be designed almost at will by changing the geometry of 
the optical devices [SSI ISO] ■ Finally, both the amplitude and the correlation length 
(down to fractions of micrometers) can be controlled accurately and calibrated using 
ultracold atoms [48] , 

Within the context of ultracold gases, important theoretical efforts have been 
devoted to disordered optical lattices which mimic the Hubbard model [HI [61] [62] 
[63] [64] . For bosons, quantum phase transitions from supcrfiuid to Bose glass and 
Mott insulator phases have been predicted [SFJ ES] and evidence of the Bose glass has 
been obtained experimentally [67] . With Fermi-Bose mixtures, the phase diagram is 
even richer and includes the formation of a Fermi-glass, a quantum percolating phase 
and a spin glass [68j ESI ED] . Effects of disorder in Bose gases at equilibrium without 
optical lattice have been addressed in connection with the behavior of the BEC phase 
transition [TT] [75] , the quantum states of Bose gases [73] [73] [75] [75] , the localization 
of Bogolyubov quasi-particles [77J [78], the dynamics of time-of- flight imaging of 
disordered BECs [79] [80], and random-field-induced order in two-component Bose 
gases [8T ] l82 ] . 

1.2. Scope and main results of the paper 

The dynamics of BECs in disordered (or quasi-disordered) potentials is also attracting 
significant attention in a quest for observing Anderson localization in non-interacting 
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BECs [55[ [53 HI or in BECs with repulsive interactions [25J [EH [El [HI] ■ Recent 
experiments have demonstrated the strong suppression of transport in expanding 
BECs in the presence of optical speckle potentials [45, 46, 48 , but this effect is not 
related to Anderson localization [15] . 

In this paper, we theoretically and numerically analyze the expansion of an 
interacting one-dimensional (ID) BEC in a disordered potential. We focus on a regime 
where the inter-atomic interactions initially exceed the kinetic energy (Thomas-Fermi 
regime), a situation that significantly differs from the textbook Anderson localization 
problem but which is relevant for almost all current experiments with disordered BECs 
[45j [46[ HZ1 SHI ESI HQ] • We distinguish two regimes that we name strong disorder and 
weak disorder respectively. 

The case of strong disorder corresponds to the situation of the experiments of 
Refs. [45;, 46,, 48] where the interaction energy in the center of the BEC remains large 
during the expansion and where the reflection coefficient from a single modulation of 
the disordered potential is of the order of unity. In this case, our numerical results 
reproduce the strong suppression of the transport of the BEC as observed in the 
experiments of Refs. [45, 46j[48]. We analyze the scenario of disorder-induced trapping 
proposed in Ref. [45] in which two regions of the BEC are identified. The first region 
corresponds to the center, where the trapping results from a competition between 
the interactions and the disorder. The second region corresponds to the tails of the 
BEC where almost free particles are multiply scattered from the modulations of the 
disordered potential. There, localization is rather due to the competition between 
the kinetic energy and the disordered potential but is ultimately due to the almost 
total classical reflection of the matterwave from a single barrier. These two effects are 
responsible for blurring Anderson localization effects [35J [35] . 

Weak disorder corresponds to a situation where the probability of large and 
wide modulations of the disordered potential is small. In this case, we show that 
Anderson localization does occur as a result of multiple quantum scattering from 
the modulations of the disordered potential. Let us briefly describe the scenario 
first proposed in Ref. 25]. Initially, the repulsive interactions are important as 
compared to the kinetic energy and to the potential energy associated to the disordered 
potential. Then, the interactions induce the expansion of the BEC and determine the 
momentum distribution of the BEC. After a time typically equal to the inverse of 
the initial trapping frequency, the interactions vanish and the momentum distribution 
reaches a steady state. Then, the BEC is a superposition of non-interacting waves 
of momentum k. Each wave localizes with its own localization length L loc (k). By 
calculating analytically the superposition of the localized waves, we show that the 
BEC can be exponentially localized or only show an algebraic decay depending on 
the correlation function of the disordered potential. In particular, due to peculiar 
long-range correlations, the BEC localizes exponentially in speckle potentials only if 
£,ia > Or, where £ in is the initial healing length of the BEC and a R is the correlation 
length of the disorder. 

1.3. Organization of the paper 

The paper is organized as follows. In section [2] we review the properties of a BEC at 
equilibrium in a combined harmonic plus disordered potential, in particular in the non- 
trivial regime where the healing length of the BEC exceeds the correlation length of the 
disordered potential. The next two sections deal with the expansion of an interacting 
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BEC in a disordered potential. Section[3]is devoted to the case of strong disorder. We 
reproduce and complete our previous results [IS] which demonstrate the suppression 
of the expansion of the BEC in a speckle potential with similar parameters as in the 
experiments of Refs. [45 , 46JHH] - The scenario of disorder-induced trapping is analyzed 
and characteristic properties of the BEC trapped by the disorder are calculated 
analytically and compared to numerical results. In particular we derive an analytic 
expression for the central density of the BEC trapped by disorder which happens to 
be characteristic of the disorder-induced trapping phenomenon and we show that the 
ultimate suppression of the expansion of the BEC is due to classical reflections from 
the large modulations of the disordered potential. We also compare these findings 
with the case of a BEC expanding in a periodic potential with similar characteristics 
as the disordered potential. Section [4] is devoted to the case of weak disorder. We 
show that Anderson localization can show up in an expanding, interacting BEC under 
appropriate conditions that are clarified. We show that the localization properties 
of the density profile crucially depend on both the momentum distribution of the 
expanding BEC and the correlation function of the disordered potential. In particular, 
in the case of a speckle potential, we find a ID effective mobility edge. We calculate 
analytically the expected localization lengths and compare our findings to the results 
of numerical calculations. Finally in section [5] we summarize our findings and discuss 
expected impacts of our work on experiments on disordered BECs. 



2. Condensates at equilibrium in a combined harmonic trap plus 
disordered potential 

2.1. Interacting Bose-Einstein condensates in a ID inhomogeneous potential 

We consider a low-temperature ID Bose gas with short-range atom-atom interactions 
gS(z) where g is the ID coupling constant. The Bose gas is assumed to be subjected to 
(i) a harmonic potential of frequency uj and (ii) an additional inhomogeneous potential 
V(z). In a finite system as considered in this work, assuming weak interactions, i.e. 
n mg/h 2 where n is the average density and m the atomic mass [HH1 IHS] ' t ne Bose 
gas will form a Bose-Einstein condensate even in low-dimensional (e.g. ID) geometries 
[59"] . Hence, we can treat the BEC within the mean-field approach [25] [35] using the 
Gross-Pitaevskii equation (CPE): 

^ + V(z) (1) 



ihdtip(z, t) 



2m 



+g\iP(z,t)\ 2 - fx 



ip(z,t), 



where fj, is the BEC chemical potential. 

In the following, we investigate the situations where the additional potential 
reads V[z) = V R v(z) with v(z) being either a disordered or a periodic function with 
vanishing average and unity standard deviation. Therefore, we have (V(z)) — and 
sj (V(z) 2 ) — (V(z)) 2 = \V R \. The sign of V R depends on the definition of the function 
v(z) and on the kind of potential one considers. For instance, in optical speckle 
potentials, the quantity v(z) + 1 is defined to be positive and V R > for blue-detuned 
laser light (case of the experiments of Refs. [45, 48, 79]) while V R < for red-detuned 
laser light (case of the experiments of Refs. [44] |46] [80] ) . For a sine-periodic potential, 
using V R < or V R > does not change the physics. See | Appendix A| for details. 



CONTENTS 



6 



2.2. The Bose-Einstein condensate wavefunction 

Here, we briefly discuss the influence of an inhomogeneous potential on the BEC at 
equilibrium in the harmonic trap. We assume that the amplitude of the disordered 
potential is smaller than the chemical potential of the BEC (V R <C /i) and that /i 3> h\o. 
This question has been investigated in details in Ref. [74]. Here, we only outline the 
results. 

At equilibrium, the BEC wavefunction is real (up to an irrelevant uniform phase) 
and is the solution of Eq. |T]) with dt4> = 0: 



fj,tp(z) 



2m 



m- (2) 



For V R — and /i ^> hoj, the kinetic term can be neglected (Thomas- Fermi regime 
[55]) and the BEC wavefunction is ipQ (z) — \J n${z) with 

[i-muj 2 z 2 /2 

n (z) = . (3) 

9 

for z such that \jl > muj 2 z 2 /2 and uq(z) = elsewhere. The density profile no(z) is 
an inverted parabola of length L TF — ^j2^,/muj 2 much larger than the healing length 
£ in = h/^/Amjl [29] (notice that this definition is different from the one of Ref. [74], 
where we used £ = \/2£in). 

In the presence of an inhomogeneous potential (Vr ^ 0), the parabolic shape of 
the density profile is perturbed. In general, the kinetic term cannot be neglected any 
longer. In particular, when £ in > a B , the short-range modulations of the potential 
V{z) induce short-range modulations of the BEC wavefunction which contribute 
significantly in Eq. ^ through the kinetic term. In order to take into account the 
effect of the inhomogeneous potential, we use a perturbative approach along the lines 
of Ref. [73]: we write i(>(z) — ^po(z) + 5i[>(z) with Sip <C tpo- The first order term of the 
perturbation series of Eq. @ is governed by the equation 

-{0 2 dl{H) + H = -^f^, (4) 

2gn 



where £° = 6n/ yl — (z/L TF ) 2 is the local healing length. Since L TF ^> (£ ln , er R ), it is 
legitimate to use the local density approximation (LDA) [2S], i-c in a region smaller 
than L TF , the quantity no can be considered as uniform. In this approximation, 
the solution of Eq. ([4]) is easily found by turning to the Fourier space. We find 
Sip(q) = -V(q)ip /2gn where 

Ha) - (5) 



and finally, 

iKz)*Mz) . (6) 

2gn (z) 

or equivalently, 

n(z) ~ n (z) - V(z)/g. (7) 
This solution justifies a posteriori the use of a perturbative approach for V R <C (J-, where 



V R = y (V(z) 2 ) — (V(z)) 2 is the standard deviation of the potential V{z). Notice that 
the equality (V(z)) — directly follows from Eq. ([5]). 
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An important consequence of the solutions © , J7J) is that the BEC wavefunction 
is only weakly perturbed by the inhomogeneous potential V(z) if V B <C A*. It follows 
from Eq. ([5]) that for £ in <C cr R , V(z) ~ V(z) and the relative inhomogeneities of the 
BEC density are 6n/n ~ V^.//i <C 1. For £ in > a R , the relative inhomogeneities are 
even smaller since all Fourier components of V are smaller than those of V. More 
precisely, the effective potential V is roughly obtained from V by suppressing the 
Fourier components with a wavelength smaller than the healing length. In other words, 
the BEC density does not follow the modulations of the bare disordered potential V(z) 
but actually follows the smoother modulations of the smoothed disordered potential 
V(z). 

Therefore, an interacting BEC at equilibrium in a disordered potential is not 
localized in the sense of Anderson. One may wonder whether this conclusion still 
holds for stronger disorder or weaker interactions, where the mcanficld approach can 
break down. This question has been addressed in Ref. 75J. It turns out that for 
very weak interactions, the Bose gas forms a so-called Lif shits glass which corresponds 
to a Fock state of various localized single-particle states. These states belong to the 
Lifshits tail of the non-interacting spectrum and are strongly trapped. Therefore, 
Anderson localization can hardly be observed unambiguously in this case. It seems 
more favorable to find evidences of Anderson localization in transport experiments of 
interacting BECs, rather than studying BECs at equilibrium in a disordered trap. 

3. Strong disorder: Suppression of the expansion of a Bose-Einstein 
condensate in a speckle potential and disorder-induced trapping scenario 

In this section, we investigate the transport properties of a coherent BEC in a 
disordered potential in the situation of the experiments of Refs. [JSJ UHl U71 05] • We 
thus assume (i) that the chemical potential of the BEC is larger than the depth of the 
disordered potential (fi > V R ) and (ii) that the correlation length of the disordered 
potential is much larger than the healing length of the BEC and much smaller than 
the (initial) size of the BEC, £ in <C (Jr <C L tf . We present numerical results which 
reproduce the suppression of the transport of the BEC in a speckle potential, observed 
in Refs. [45, 46, 48J, and discuss a scenario to explain this phenomenon. In addition, 
we compare the observed behavior to the case of a periodic potential with similar 
characteristics. 

3.1. Expansion of an interacting BEC in a speckle potential 

In order to induce transport, we start from a BEC at equilibrium in the harmonic and 
disordered potentials (see Sec. H]). At time t = 0, we suddenly switch off the trapping 
harmonic potential, keeping the disordered potential. This process is similar to the 
one used in Refs [45], |46l |48l [85] . The evolution of the BEC is thus governed by the 
CPE (fl]) with iij = and the initial condition corresponds to the TF wavefunction 
discussed in Sec. 12.21 

The time evolution of the root mean square (rms) size of the BEC, Az(t) = 
yj (z 2 ) — (z) 2 , as obtained from the numerical integration of the time-dependent 
CPE |T]), is plotted in Fig. [T] for several amplitudes V R of the disordered potential. 
In the absence of disorder, the interacting BEC expands self-similarly as predicted by 
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Figure 1. (color online) Time-evolution of the rms-size of the BEC wavefunction 
evolving in the disordered potential V for several values of the amplitude Vr. 
The (black) dashed line is the theoretical prediction of the scaling theory dlOI I 
with a vanishing disordered potential. Here, we have used Ur. = 0.012L T f and 
&„ = 5.7 x 10- 4 L TF . 



the scaling approach [501 HI] : 

,, s ib[z/b(t),0] ( mz 2 b(t)\ 

tp(z,t) ~ yv '111 ' exp ^-ttV (8 
y/btf) V 2hb ^ J 

where b(t) is the scaling parameter which is governed by the equation 

8(f) = tu 2 /b 2 (t) (9) 

with the initial conditions b(t = 0) = 1 and b(t = 0) = 0. Integrating these equations, 
we find 

VHt)(b(t) - 1) + H-v/ftC*) + VKt) - 1] = V2ut, (10) 

which asymptotically reduces to a linear expansion at large time, b(t) ~ \/2wt. The 
numerical calculations agree with this expression as shown in Fig. [TJ 

The situation is significantly different in the presence of disorder. In this case, the 
initial BEC wavefunction is the usual Thomas-Fermi inverted parabola perturbed by 
the disordered potential T4J. For t < l/u>, the scaling form ((8|) is still a good solution of 
the GPE and, according to the scaling theory [9Ql [91] the BEC wavefunction expands. 
For larger times and small amplitudes of the disordered potential (V R < 0.1/i), the 
effect of disorder on the expansion observed in the numerical calculations is small and 
the BEC expands by about one order of magnitude for lot =10. For larger amplitudes 
of the disorder (V R > 0.15/x), the expansion of the BEC stops after an initial expansion 
stage described above. This effect signals the localization of the BEC wavefunction 
due to the presence of disorder. 

Important information can be obtained from density profiles of the localized BEC. 
For instance, density profiles corresponding to a single evolution are plotted at two 
different times in Fig. [21 From these, it appears that the localized BEC is made of 
two distinct parts: a static dense core and fluctuating dilute tails (see also Fig[3]). In 
particular, the small fluctuations of Az observed in Fig. Q] are due to the contribution 
of the tails of the BEC that still evolve while the core of the wavefunction is localized 
(see section [3T2T) . 
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Figure 2. (color online) Density profiles of the BEC for Vr = 0.2/^, or = 
0.012L T f and £ in = 5.7 X 10 _4 Ltf (solid red lines) for two different values of 
the expansion time r and expected Thomas-Fermi profiles in the absence of a 
disordered potential (dashed green lines). We also show the disordered potential 
normalized so as to be homogeneous to a density (V/g; dotted blue line). 



It is worth noticing that the BEC expansion stops for amplitudes of disorder 
significantly smaller than the typical energy per particle in the initial BEC: V R < ji. 
This suppression of transport is phenomenologically similar to what is expected from 
Anderson localization [8l [9l [10]. Strictly speaking, Anderson localization relies on 
the existence of localized single-particle eigenstates and on the subsequent absence of 
diffusion [8]. However, we have stressed that the presence of predominant inter-atomic 
interactions dramatically changes the picture [45] . On one hand, repulsive interactions 
are expected to reduce the localization effect (74J [92] . During the initial expansion 
of the BEC, the interaction energy greatly dominates over the kinetic energy in the 
center of the BEC so that no Anderson-like localization effect is expected in this region. 
On the other hand, although the particles in the tails are weakly interacting due to 
the small density, the initial interactions determine their typical energy as the initial 
expansion stage converts the interaction energy into kinetic energy. We will see that, 
for strong disorder as considered here, the modulations of the disordered potential will 
ultimately stop the expansion of the dilute tails, masking any Anderson-like effect (in 
the case of weak disorder however, Anderson localization can be obtained in this region 
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as discussed in Sec. HJ). In the following, we detail the scenario of disorder-induced 
trapping outlined above and first proposed in Ref. [45j . 



3.2. Scenario of disorder-induced trapping 

The dynamics of the BEC in the disordered potential is governed by three different 
forms of energy: (i) the potential energy associated to the disordered potential, (ii) 
the interaction energy and (iii) the kinetic energy. It is thus useful to evaluate and 
compare the kinetic and interaction energies to understand the behavior of the BEC 
in the disordered potential. To this end, notice first that it follows from the initial 
expansion of the BEC that the fast atoms populate the tails of the expanding BEC 
while the slow atoms stay close to the center. In addition, notice that, except for very 
small amplitudes of the disordered potential and subsequent long expansion times, the 
density in the core of the BEC remains large whereas it drops to zero in the tails (see 
Fig. [2]). We thus distinguish two different regions of the BEC: (i) the core where the 
density is large and the interaction energy is dominant and (ii) the tails where the 
density is small and the kinetic energy dominates. The behavior of the BEC turns out 
to be completely different in these two regions [35] . 



3.2.1. Quasi-static Thomas- Fermi profile in the core of the BEC - For the sake of 
clarity, we define the core of the BEC as half the total size of the initial condensate: 
— L TF /2 < z < L TF /2 and call 

».(*) = 7^- [ +LTF 2 dz\i>(z,t)\ 2 , (11) 

L TF J-Ltf/2 

the average BEC density in the center. In particular, at time t = 0, due to the 
parabolic envelope resulting from the harmonic trap, we find 

11m , 
12 g 

in the absence of disorder but also in the presence of a self-averaging disordered 
potential]]. 

During the initial expansion stage, the average density in the core, n c , slowly 
decreases and the parabolic envelope disappears. Since the interaction energy 
significantly exceeds the kinetic energy, we expect the local density \ip(z,t)\ 2 to follow 
almost adiabatically the instantaneous value of n c (t) approximately in the Thomas- 
Fermi regime so that 

\^(z,t)\ 2 ^n c (t)-V(z)/g. (13) 

In order to check this prediction, we plot in Fig. [3^ the result of the numerical 
integration of the CPE (JTJ) for the density profile in the central region of the BEC 
during the evolution in the disordered potential at two different times, together with 
a plot of the analytical expression (|13p . In particular, two properties are of special 
interest here. First, we observe that the time-dependent fluctuations of the density 
profile are significantly smaller than the modulations of the disordered potential 
V(z)/g. Second, the density profiles are in good agreement with Eq. (fT3|) . This 

| In the context of disordered systems, a quantity is said to be 'self-averaging' when it verifies 
the principle of 'spatial ergodicity'. In other words, it means that the average over realizations of 
the disordered potentials of a relevant quantity, F, equals the corresponding spatial average: i.e. 
(F) ~ i f Q L dzF(z). 
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Figure 3. (color online) Density profiles of the BEC during the evolution in the 
disordered potential at different times for Vr = 0.2^ in the core (a) and in the 
tails (b) of the BEC. Both are magnifications of the plots of Fig. [2] The solid 
(red online) and dashed (blue online) lines correspond respectively to the times 
cut = 10 and ut = 20 of the same evolution and the dotted (purple online) line 
corresponds to Eq. I|13j l with n c as a fitting parameter. Notice the different scales 
in the two figures. 

observation supports the scenario of an adiabatic decrease of the density in the center 
of the BEC. The value of n c at the end of the expansion turns out to be characteristic 
of this scenario. In the following, we show that n c can indeed be computed from the 
statistical properties of the disordered potential. 

The expansion of the core of the BEC in the disordered potential stops when 
the condensate fragments, i.e. when the effective chemical potential in the center of 
the BEC (~p, = gn c ) decreases down to the value of typically two large modulations 
of the disordered potential. At this time, the energy per particle in the core of the 
BEC becomes too small to over-pass the potential barriers and the core of the BEC 
gets trapped between these large modulations. This scenario allows us to determine 
the final value of the average density n c in the core of the BEC. Let us call N p ^ kB (V) 
the number of maxima of the disordered potential in the central part of the BEC 
(— i TF /2 < z < L TF /2) with an amplitude larger than a given value V and assume 
that it can be computed from the statistical properties of the disordered potential. 
The density in the center of the BEC after the trapping has occurred thus corresponds 
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Figure 4. (color online) Average density n c in the core of the BEC trapped 
by the disorder versus the amplitude of the disordered potential Vr for different 
values of the correlation length cr R and comparison to Eq. IU4I , The horizontal 
(red online) line corresponds to the saturation limit n c = llfj,/12g. 



to the maximum value of n c below which two modulations of V in average are present 
in the center of the BEC. This is simply computed by solving for N pcB . ks (V = n c g) = 2. 
Although, this scheme is general, it appears clearer when applied to a case where 
^Vpeaks(V') can be explicitly computed. Let us now consider the case of a speckle 
potential [SSl[5D] with V R > 0. It is shown in |Appcndix A| [see Eq. (|A.5[) ] that in this 



exp 



-0 



where a ~ 0.30 and /3 ~ 0.75. From this, we 



case N P ^ S {V) ~ a (^) 

easily find that the final density of the core of the BEC is n ~ ^ (~g) ^ n 
addition, we notice that the final density cannot exceed the initial density as resulting 

fi/g. In 



2<T R 



In 



fro 



In 



2cr R 



m an expansion. Therefore Eq. (fll| is valid only for ^ (~g) 
the opposite situation, the BEC is already multiply fragmented at t — and the final 
density saturates at n c ~ [see Eq. (fT2j) ]. In summary, we expect that the average 
density of the BEC trapped by the disorder is 



1 

m 



in 



aL-r 



2a R 



llfi 
125 



(14) 



In order to check Eq. (Ti"4"|) , we have extracted the averaged central density [see 
Eq. (flT|) ] from the wavefunctions ip calculated numerically for several amplitudes V R 
and correlation lengths <r B of the disordered potential. In Fig. [U we plot n c as a 
function of V R for several a R together with the prediction (|14j) . The results show that 
Eq. ([T4^) provides a good estimate of the final density n c in the core of the BEC. 
In particular, for small amplitudes of the disorder, n c grows linearly with V R with a 
coefficient in agreement with Eq. ([14]) up to about 10 %. For larger amplitudes of the 
disorder, n c saturates below \\\ijYlg as expected. 

This behavior agrees with experimental results for a blue-detuned speckle 
potential (V R > 0) 48J . It is worth noticing that our scenario is expected to apply also 
to the case of a red-detuned speckle potential (V R < 0) as used in Ref. [IB]. In this 
case, the fragmentation occurs when fx = \V R \ independently of the correlation length 
of disorder (if £ in <C cr n ). Then, the fragmented BEC is trapped in the small wells of 
the disordered potential with a typical size a R and with a central density n c ~ |Vr|/<? 
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(independent of cr R ). Instead, for a blue-detuned speckle potential as investigated 
above, the BEC is trapped between large modulations that may be separated by a 
distance much larger than o~ B . As a consequence the final density at the center is 
expected to be significantly larger. This is confirmed by the numerical results of 
Ref. [93]. 

3.2.2. Strong reflections in the tails of the BEC - The situation is completely 
different in the tails of the BEC. Due to the small atomic density, the kinetic energy 
now dominates over the interaction energy. The tails are populated by fast moving, 
weakly interacting atoms that undergo multiple scattering from the modulations of 
the disordered potential. Ultimately, the trapping of these atoms results from almost 
total classical reflection on a single large modulation of the disordered potential with 
an amplitude exceeding the typical energy of a single particle [45j [46] . This scenario 
is supported by the density profiles plotted in Fig. [2] where one can observe a sharp 
drop of the atomic density at the edges of the BEC {i.e. at positions z m j n — — 7L TF 
and Zmax — 8L TF in Fig. \2§. Notice that significant drops correspond either (i) 
to modulations of the disordered potential larger than the initial chemical potential 
/i (e.g. at z mm ~ — 7L TF ) or (ii) to a concentration of weaker barriers (e.g. at 
^min ~~ 3.5-L TF ). 

In contrast with the situation in the core of the BEC, it is expected that (i) the 
density profile does not show a Thomas-Fermi shape and (ii) the local density is not 
stationary. Both properties agree with our numerical results as shown in Fig.[3jD where 
we plot a magnification of a small region corresponding to the tails of the BEC of Fig. [2] 
In particular, the shorter modulations of the wavefunction observed in Fig.[3]3 are due 
to the kinetic energy of the particles in the tails. This statement is corroborated by 
the calculation of the energy per particle, e. Due to energy conservation, the energy 
can be computed at the initial time t = (i.e. right after releasing the BEC from the 
trapping potential), e = J dz \ip(z)\ 2 [V(z) + g\ip(z)\ 2 /2] . Using Eq. ©, we easily 
find that 

(15) 



2//. 

T 



8 U 



The disordered potential perturbs the energy per particle only at second order in V K /fi, 
and, for V B <C /i, we have e oc fi. From this, we expect that the typical wavelength A 
of the fluctuations in the tails would be of the order of the healing length in the initial 
condensate, so that A/2n ~ £ in . This is confirmed by the properties of the momentum 
distribution of the BEC which show two sharp peaks located around p ~ ±h/£ in . 

3. 3. Expansion of a condensate in a periodic potential 

Up to this point, our analysis has shown how the competition between inter- atomic 
interactions and disorder (in the center of the BEC) or between kinetic energy and 
disorder (in the tails of the BEC) can strongly suppress the coherent transport of an 
interacting matterwave in a disordered potential. A natural extension of our analysis 
is to compare these findings to the situation in a periodic potential with similar 
characteristics (see | Appendix A I. In the case of a periodic potential, no suppression of 



transport is expected as no large peak can provide a sharp stopping of the expansion, 
and obviously, no 'a la Anderson' localization should occur. 
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Figure 5. (color online) a) Time-evolution of the rms-size of the BEC 
wavefunction evolving in a periodic potential for several amplitudes Vr. and for 
A = O.IILtf- The theoretical prediction corresponding to Eq. (1101 in free space 
is also shown (black dotted line). The inset shows the velocity of the expansion 
of the BEC together with the theoretical estimate (1171 . b) Density profiles of the 
BEC in the harmonic trap and after an expansion time in the periodic potential 
of t = 40/u; for Vr = 0.2/x. 



Numerical results for the expansion of the BEC in the periodic potential described 
in |Appcndix A.3| are shown in Fig. [SK The difference with the case of a disordered 
potential (see Fig. [1} is striking: as expected, the BEC now expands linearly with 
time with an asymptotic expansion rate that decreases when the amplitude V R of the 
periodic potential increases. A detailed analysis shows however that the transport of 
the BEC in a periodic potential and in a disordered potential share some properties 
for the parameters used in this section, as we discuss below. 

Again, important information is contained in the density profiles such as the ones 
plotted in Fig. [5b (to be compared to Fig. [3] which corresponds to the disordered 
case). Initially (cut — 0), the density profile follows the modulations of the periodic 
potential modulated by the parabolic envelope associated to the harmonic trapping 
[74] , During the initial expansion stage, the density in the center decreases slowly and 
follows adiabatically a Thomas-Fermi shape with a slowly decreasing instantaneous 
chemical potential ~p. Then, the evolution of the center stops when the chemical 
potential JL exactly matches the potential depth, i.e. when the BEC fragments. This 
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Figure 6. (color online) Average density in the center of the BEC trapped in 
the periodic potential versus the lattice depth Vr and for two lattice spacings A. 
The theoretical prediction H16H is also plotted (green dashed line). 



is similar to the disordered case. However, in the case of a periodic potential, it is a 
deterministic process which appears when 

As shown in Fig. [6l this formula provides a very good value of the average density in 
the center of the BEC trapped in the periodic potential. Remarkably, Eq. (fT6|) shows 
that n c does not depend on the lattice spacing A as confirmed by the numerical results 
shown in Fig. [S] This is different from the case of a blue-detuned speckle potential 
(Vr > 0) where n c has been shown to depend explicitly on the correlation length of 
the disordered potential [see Eq. (|14p and Fig. 0] . 

During the subsequent evolution, a part of the BEC is thus trapped in the 
center while the tails still expand as shown in Fig. [5J Let us focus now onto the 
tails of the BEC. To do so, let us first write the BEC wavefunction ifj = tpc + i/Jw 
where 4>c and tpw account for the center (|z| < L TF ) and for the tails (\z\ > L TF ) 
respectively. As the supports of ipc and -0 W are spatially separated, we have 
Az 2 = J dz z 2 1 -0c | 2 + / dz z 2 |-0w| 2 - The center gets trapped after a transient time 
so that / dz z 2 \ipc\ 2 tends to a constant, Azq, at large times. In contrast, the tails 
expand so that their density decrease. After a typical time 1/to, a substantial part of 
the interaction energy is converted into kinetic energy and the interaction term can be 
neglected for the subsequent dynamics. Wc also neglect the periodic potential which 
has a small amplitude compared to the typical energy per particle in the tails. Now, 
in free space, Az 2 ^ = J dz z 2 \if; w \ 2 ~ i^"^ 2 at large times where E w is the total 
(kinetic) energy in the tails of the BEC. Due to the conservation of the total energy 
during the expansion, we have Ne = Eq + E w where e is the energy per particle 
given by Eq. (fT5|) . and the energy in the center of the BEC, Eq — V 2 L TF /g, is easily 
computed from the Thomas- Fermi profile in the center of the BEC [see Eq. (fT3|) ] . We 

2"| 



finally find that Ew /N ■ 



2yi 

5 



1 



15 f EelY 
4 U J 



so that Az 2 ~ Azq +v 2 t 2 at times larger 
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than 1/lo with 





) 



21 V2 



— Wi. 



TF 



(17) 



In the absence of disorder (V R = 0), Eq. (fTT|) is consistent with the scaling theory 
[Eq. (10])] for which Az(t) = b(t)L TF /VE. In the presence of disorder, it provides a 
reasonable agreement with the numerical findings as shown in the Inset of Fig. [5^- We 
attribute the discrepancy at the largest values of V K to the main two approximations 
that have been used. First, a strict separation between the tails and the center 
at z = L TF has been used to compute E-w and we have thus neglected the small 
intermediate region. Second, the interaction of the atoms with the periodic potential 
is expected to increase the inertia of the expanding gas and this should contribute to 
slightly lower the expansion velocity compared to the prediction (|17|) . 

4. Weak disorder: Onset of Anderson localization in the expansion of a 
condensate 

In section 13] we have shown that in the experimental conditions of Refs. 45, 46j 08] , 
Anderson localization effects are blurred in an expanding, interacting BEC owing (i) to 
important repulsive interactions in the center of the BEC, and (ii) to strong reflections 
from single barriers of the disordered potential in the tails of the BEC. Both effects 
are related to the presence of large modulations of the disordered potential. It thus 
appears necessary to work in a parameter range where the probability of single large 
modulations of the disorder is negligible in order to observe unambiguous Anderson 
localization of an expanding BEC. 

In this section, we work within the regime of weak disorder [a precise definition 
is given below, see Eq. (|2T|) ]. Following the theory of Ref. [25], we show that in 
this situation, Anderson localization of an interacting BEC can be observed under 
appropriate conditions that we identify precisely [25] . We consider both cases of an 
impurity model of disorder and of a speckle potential. In particular, for a speckle 
potential, we show that the long-range correlations induce a ID effective mobility 
edge, i.e. strong exponential localization is obtained only for £ in > o~ B . 

4-.1. General model of Anderson localization of an expanding Bose-Einstein 



Expansion of a Bose-Einstein condensate - Let us examine again the expansion of 
the BEC in the disordered potential (see section [XT]) . For weak disorder, the initial 
interaction energy strongly exceeds the potential energy associated with the disorder 
so that the first stage of expansion of the BEC is hardly affected by the disorder. For 
instance, the numerical results of Fig. [1] for V n = 0.5/i and V R — O.l/i confirm this 
assertion for durations of expansion up to about t ~ 10/w. Within this time window, 
the momentum distribution of the expanding BEC can thus be approximated to that 
of a BEC expanding in free space [see Eq. |J5])]. Calculating the Fourier transform 
of the scaling solution for interacting BECs expanding in free space, ip(z,t), using 
the stationary phase approximation (valid for t ^> h/ff), we find the momentum 



condensate 



CONTENTS 



17 



distribution 



V{k,t) 



37V£ ln /4 

Vi - VW) 



x e 



where 8 is the Heaviside step function. Since b(t) 
reaches a steady-state at times t 3> l/u>: 

3iV£ tI 




(18) 



\/2u)t, the momentum distribution 



V{k) 



x [i - (ku 2 ] x e [i - (fc^n) 2 



(19) 



An important feature of the momentum distribution (|19p is that it has a high- 
momentum cut-off at k c = 1/Cin (see Fig. [7^,). 

For t l/u), almost all the initial interaction energy is converted into kinetic 
energy. Neglecting the effect of disorder at this stage, this property can be obtained 
from the scaling solution ([5]) [SU1 US]- We find that the interaction energy is 
E iat (t) ~ E int (0)/b(t). Then using the property of conservation of the total energy 
during the expansion, we find that the ratio of the kinetic energy to the interaction 
energy is E kin (t) / E iat (t) ~ b(t) — 1 which is much larger than unity for t ^> l/u. It 
follows from this analysis that for times typically larger than l/u>, the expanding BEC 
is a coherent superposition of almost non-interacting plane waves of momentum k: 

F dk ~ 

iK*,f)= / -=v(MK fez , (20) 



2tt 

the momentum distribution T>(k) — \ip(k,t)\ 2 being stationary and determined by the 
interaction-driven first expansion stage [25) . 



Anderson localization of quantum single particles in a correlated disordered potential - 
Therefore, the interaction of each fc-wave with the disordered potential can be 
treated independently. According to the Anderson theory [8], the k- waves will 
exponentially localize as a result of multiple scattering from the modulations of the 
disordered potential. In other words, each component e lkz in Eq. ([2H| will become 
a localized function 4>k(z), characterized by an exponential decay at large distances: 
ln|(/)fc(z)| ~ — 7(fc)|z|, where j(k) — l/L loc (k) is the so-called Lyapunov exponent, 
and L loc (k) is the localization length. The Lyapunov exponent can be calculated 
analytically in a correlated disordered potential using the phase formalism approach 
[M] (see also | Appendix B[ ). At the lowest order of the Born expansion, which is valid 
provided that j(k) <C k, i.e. for 

V R a R « — (fc^.) 1 / 2 , (21) 

TO 

we find the Lyapunov exponent 

7 (fc) ~ ^ (fca R ) 2 c(2fc CTR ), (22) 

where E = h 2 k 2 /2m, and c(n ) = J ^=c{u )e lKU is the Fourier transform of the reduced 
correlation function c{u) (see |Appcndix A\ . A plot of the Lyapunov exponents versus 
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Figure 7. (color online) a) Stationary momentum distribution of an expanding, 
interacting Bose-Einstein condensate at large times [t > X/u>). b) Lyapunov 
exponent of single-particles of energy E as a function of the momentum k = 
V2mE/h. 



the momentum k is shown in Fig.[7j3 for a speckle potential and for a Gaussian impurity 
model. 

Deviations from a pure exponential decay of <pk turn out to be important here. 
They can be obtained using diagrammatic methods 95, 96], and one finds an integral 
formula for the average density of each fc-wave: 

{\4>k{z)\ ) = I duu smh(7TL^ ' 



2 J \ 1 + cosh(7ru) t 

x exp{-2(l + w 2 ) 7 (fc)|z|}, (23) 



where j(k) is given by Eq. ([22]) . Notice that at large distances (\z\ ^> l/j(k)), Eq. (|23|) 
reduces to 

M .f^=)x^Ei^M, (24) 



64 



3/2 



so that the exponential decay of the density of the localized single-particle states is 
corrected by an algebraic decay l/|z| 3 / 2 . 



Anderson localization of the Bose-Einstein condensate - In the regime of weak 
disorder defined by condition (|21|) . the Anderson localization transforms each plane 
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wave e lkz which appears in the superposition (|20p into the localized wave <p z (z) — 
r(z) sin[0(z)] where 9(z) ~ kz and r(z) is a slowly decaying envelope (see |Appcndix B[ ). 
Therefore, the Fourier transform of 4>k{z) is peaked around k. It follows that 
the interaction of the fc-wave with the disordered potential only weakly affects the 
momentum distribution of the BEC. Hence, once each independent fc-wave is localized, 
the density of the BEC is given by the equation 

2\ 



n (z) 



dk » 

ip(k,t)<f> k (z) 



2tt 



(25) 



with \ip(k,t)\ 2 ~ T)(k). Assuming that the phases of the functions 4>k{z) 1 which are 
determined by the local properties of the disordered potential and by the evolution 
time, are random, uncorrelated functions for different momenta, i.e. (<f>p(z)(pk(z)) ~ 
27T(5(fc — fc') , the density of the BEC reduces to 

/•OO 

n (z)~2 dkV{k){\<j> k (z)\ 2 ) (26) 
Jo 

where we have used the properties T)(k) = T>(—k) and (\4>k{z)\ 2 ) — (|0_fc(z)| 2 ). This 
formula is transparent and contains the main ingredients of the Anderson localization 
of an interacting BEC in a disordered potential. In a first stage, the interactions drive 
the expansion of the BEC and determine the momentum distribution T){k). In a second 
stage, the interactions vanish and the BEC is formed of a superposition of plane waves 
of energies E = H k /2m. Then, each fc-wave localizes with its own localization length 
Li oe (k) — l/j(k). We will show in the next paragraphs that the precise localization 
properties of the BEC which are determined by the integral (|26|) , strongly depend on 
the correlation function of the disordered potential. This is reminiscent of the strong 
dependence of the single-particle Lyapunov exponent 7(fc) on the correlation function. 



4-2. Anderson localization of an expanding Bose-Einstein condensate in an impurity 
model of disorder 

Let us consider the case of the impurity model of disorder described in the 
|Appcndix A.2| It is made of a series of Gaussian peaks of width w and amplitude 
Vb, all identical and spread randomly along the z axis. This potential is generic and 
in the limit w — > with Vqw fixed, we recover the widely used 5-correlated (i.e. un- 
correlated) disorder made of a random series of <5 peaks used in a number of theoretical 
investigation of disordered systems [94] . This Gaussian impurity model of disorder can 
also be implemented using the so-called impurity atom technique with ultracold atomic 
gases [Sni[57]. In this case, the Fourier transform of the reduced correlation function 
reads 

c(k) = V^V2exp(-K 2 /4), (27) 

and the amplitude and correlation length are V R = WWVq and cr R = 2w respectively 
(see |Appendix A.2[ ). Inserting Eq. (J27J) into Eq. (|2"2"]l. we find 

7(fc)- ^!f 2 gR exp[-(fca R ) 2 ] (28) 

which is plotted in Fig. \7jp (blue, dotted line). 

Using Eqs. (fT§)) . (f2"5)) . ([2l)j) . (|2"gf . we now calculate the density profile of the localized 
BEC. Since the density profile no(^) is the sum over fc of the functions (\(f>k(z)\ 2 ) 
which decay exponentially with a rate 27 (fc), the long-tail behavior of no(z) is mainly 
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determined by the fc-components with the smallest j(k), i.e. those with k close to 
the high- momentum cut-off fc c = l/£ in . Therefore, integrating in Eq. (J^H) we limit 
ourselves to the leading order terms in Taylor series for T>(k) and j(k) at k close to 
fc c . We find 

exp{-27 cff |z|} 

n °( z ) K up ( 29 ) 

where 7 eH = 7(fc = l/£ in ). (30) 

This means that the Anderson localization of an expanding, interacting BEC 
occurs, provided that the disordered potential is weak enough. In the case of a 
Gaussian impurity model of disorder, the density profile shows an exponential decay 
with the effective Lyapunov exponent equal to the one of a single particle of momentum 
k = k c = l/f Jn [see Eq. (J3DJ)], i.e. 

7 rff = ^(^) 2 K/6Jcxp[-(a R /U 2 ]. (31) 

This is a clear characteristics of Anderson localization of the BEC which can be 
observed in experiments on ultracold atoms using direct imaging techniques. 

4-3. Anderson localization of an expanding Bose-Einstein condensate in a speckle 
potential 



Let us examine now the case of a speckle potential (see Appendix A.l[ ) for which, in 



ID, the Fourier transform of the reduced correlation (|A.4|) function reads 

c(k) = yf^~/2{l - k/2)6(1 - k/2), (32) 

where is the Heaviside step function. The case of a speckle potential is particularly 
interesting for two reasons. First, it corresponds to the model of disorder used in 
almost all present experiments with disordered BECs [UJ HFJ SSI EH HH1 HO] ■ Second, 
we will see that speckle potentials offer much richer situations than the impurity model 
discussed above due to peculiar long-range correlations [25] . 

One important feature of the speckle potential is the fact that the Fourier 
transform (|32| of the correlation function has a finite support. After Eq. (f22|) . it 
results that the Lyapunov exponent vanishes for k > l/a a , i.e. that strong Anderson 
localization occurs for non-interacting fc-waves only for k < l/cr R [25] - In other words, 
there is a ID mobility edge at l/er R in the Born approximation. Strictly speaking, 
higher orders in the Born expansion may provide a non- vanishing Lyapunov exponent 
for k > 1/ctr. However, we have shown using direct numerical calculations that the 
localization length (inverse Lyapunov exponent) for k > l/cr R strongly exceeds typical 
sizes of ultracold atomic samples, so that we can consider k = 1/<t r as an effective 
mobility edge in our problem [25j . 

It follows that a part of the expanding BEC (i.e. its Fourier components with 
k > l/c R ) expands to infinity while all the Fourier components with k < l/cr R localize 
exponentially with the fc-dependent Lyapunov exponent 

7TTT1 ^ fT 

7(*0 = {l " k<jR)Q{1 " kanl (33) 

found by inserting Eq. $S2§ into Eq. f2"2"|). Equation ([33} is plotted in Fig. [7b (solid, 
red line) and show that for a ID speckle potential, the high-momentum cut-off 
fc c = min{l/^ in , 1/ctr.} in the integral formula ([26)) for the BEC density is twofold. 
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The cut-off k = 1 /£ in is related to the momentum distribution of the expanding BEC 
and is due to the initial atom-atom interactions, while the cut-off k = 1/<t r is related 
to the correlation function of the ID speckle potential and is due to the peculiar finite 
range correlations of the disordered potential. Now, two very different situations must 
be distinguished [25] . 

For £ in > <7 R , the high- momentum cut-off k c is provided by the momentum 
distribution. In this case all non-interacting functions (\(f>k(z)\ 2 ) are exponentially 
localized with a finite Lyapunov exponent, "f(k) > 7(l/£ in ) > 0. This situation is then 
similar to the case of the Gaussian impurity model and, integrating Eq. (|26[) . we find 

exp{-27 cff |z|} 

n °( z ) — — ( 34 ) 

where % a = -f(k = l/£ in ). (35) 

Finally, the BEC density profile is exponentially localized with the effective Lyapunov 
exponent 

7e« = ^ (jf\ K/e in )(l - - Or/^)- (36) 

For £ in < er R , the situation is completely different. In this case, the cut-off k c 
is provided by the correlation function, and since j(k = l/c R ) = 0, the relevant 
Lyapunov exponents (^f(k) for all k < 1/<7 R ) do not have a finite lower bound. Then, 
integrating Eq. (f2"u]l . we find that the BEC density profile is not exponentially localized 
but rather shows an algebraic decay |25j : 

n (z) oc -L. (37) 
\A 

We now present numerical results performed within the Gross-Pitaevskii 
approximation for the expansion of a BEC in a speckle potential [25] . The inset 
of Fig. [5] (upper panel) shows that the expansion is strongly suppressed and for long 
times, the BEC density profile is localized as shown in Fig. [5] (upper panel). Let us 
discuss now the behavior of the tails. 

For £ in > er R , the density profile obtained numerically is found to be exponentially 
localized. In addition, fitting the function Uq{z) oc exp{— 27 cff |z|}/|,2| 7 / 2 to the 
numerical results with the amplitude and 7 cff as fitting parameters, we find that the 
results for j en are in excellent agreement with the prediction (I36|) as shown in Fig. [8^ 
(lower panel). 

For £ in < u R , we find that the density profile decays algebraically. We fit the 
function no(z) oc l/|z|' 3c/ - f with the amplitude and (3 e ff as fitting parameters, and we 
find that /3 e tf ~ 2 in agreement with the prediction (|37|) as shown in Fig. [HJa (lower 
panel). 



5. Conclusion and perspectives 

In summary, we have theoretically investigated the localization of an expanding ID 
BEC with repulsive atom-atom interactions characterized by the initial healing length 
£ in in a disordered potential of finite correlation length er R . We have restricted our 
study to the regime where the initial interactions of the trapped BEC dominate 
over the kinetic energy and the disorder, a situation relevant to almost all current 
experiments with disordered BECs [H H3 HO |H [751 ED]- When the BEC is 
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Figure 8. (color online) Upper panel: Density profile of the localized BEC in a 
speckle potential at t = 150 /ui. Shown are the numerical data (black points), the 
fit of the result from Eqs. H19I I, H23H and J26D [red solid line], and the fit of the 
asymptotic formulas (134 6 . J 3 5 1 ) [blue dotted line]. Inset: Time evolution of the rms 
size of the BEC. The parameters are Vr = 0.1/j, £i„ = O.OILtf, and cr R = 0.78£i„. 
Lower panel: a) Lyapunov exponent 7 e ff in units of 1/Ltf for the localized BEC 
in a speckle potential, in the regime £ in > or. The solid line is 7(l/£i n ) from 
Eq. (1 36 D - b) Exponent of the power-law decay of the localized BEC in the regime 
£in < or. The parameters are indicated in the figure. 



released from the trapping potential while keeping the disordered potential on, we 
find a strong suppression of the expansion, similar to earlier experimental observations 
[I5I SHI HE] ■ We have shown that this localization effect has completely different causes 
depending on the strength of the disorder. 

Strong disorder - The case of strong disorder corresponds to the situation where 
several modulations of the disordered potential individually are strong-enough to 
induce an almost total reflection of noninteracting particles of energy equal to the 
typical expansion energy per particle of the BEC. In particular this case is relevant 
to the experiments of Refs. [HI 0S1 HE] and to the numerics of Refs. g5[ [M] W\ ■ In 
this case, the localization results from a disorder-induced trapping [45] , whose scenario 
involves two processes: (i) the fragmentation of the core of the BEC on one hand, and 
(ii) classical total reflections from single large modulations of the disordered potential 
on the other hand. In the core of the BEC, the interactions remain important during 
the initial expansion stage. The BEC is in a quasi-static Thomas-Fermi regime with an 
effective chemical potential which first slowly decreases during the expansion. When 
the BEC fragments, the expansion of the core of the BEC stops. In the tails of the 
BEC, the interactions are negligible and the particles undergo multiple scattering from 
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the modulations of the disordered potential but the expansion is ultimately stopped 
by single large modulations. Hence, in the case of strong disorder, the localization is 
not related to Anderson localization. 

Weak disorder - In the case of weak-enough disorder, the probability of modulations 
of the disordered potential such that the reflection of a particle of typical energy /i 
approaches unity is negligible, and Anderson localization can show up in an expanding 
BEC. The scenario is then as follows [35] • In a first stage, the interaction energy 
dominates over both the kinetic energy and the disorder, and drives the initial 
expansion of the BEC. After a typical time of 1 /uj where to is the frequency of the initial 
trapping potential, the interaction energy vanishes and the momentum distribution 
of the expanding BEC becomes stationary. At this stage, the interactions can be 
neglected and the BEC wavefunction is a superposition of (almost) non-interacting 
waves of momentum k. Each fc-wave Anderson localizes with its own fc-dependent 
localization length L lac (k). The density profile of the localized BEC is thus the 
superposition of the localized fc-waves and strongly depends on the correlation function 
of the disordered potential. 

The case of speckle potentials is particularly interesting as the Fourier transform 
of their correlation function has a finite support. It follows that the localization of 
the expanding BEC is exponential only for £ in > er R (in the lowest order of the Born 
expansion, see Sec. 2]). In the opposite situation (£, n < er R ), the density profile of 
the BEC decays algebraically as l/|z| 2 . Therefore, for speckle potentials, there is 
an effective mobility edge at £ in = er R for the Anderson localization of an expanding, 
interacting BEC [25] . 

Perspectives - Our results suggest that the ID Anderson localization can be observed 
in an interacting BEC (initially in the Thomas-Fermi regime) expanding in a 
disordered potential in experiments similar to those reported in Refs. J35J SSI HE]- 
We stress that special attention should be paid to using weak-enough disorder to 
allow the dilution of the BEC during the first expansion stage and to avoid strong 
reflections from large modulations of the disordered potential. In addition, we have 
shown that the correlation function of the disordered potential plays a crucial role for 
the localization properties of the BEC. 

However, a couple of challenges have to be taken up to observe Anderson 
localization in current experiments with expanding BECs. Indeed, in addition to 
being able to produce long-enough expansions and to measure very small densities, it 
appears that both disorder and interactions have to be carefully controlled. 

In this respect, using disordered potentials created by optical speckle patterns is 
particularly promissing. On one hand, from a practical point of view, the correlation 
functions of speckle potentials are very well controlled and can be designed almost 
at will. This allows for a direct comparison between experimental observations and 
theoretical predictions. On the other hand, for a ID speckle potential, the Fourier 
transform of the correlation function has a cut-off at k = 1/2<t r which induces an 
effective mobility edge at k — 1/<t r for single-particles and, correspondingly, at £ in = cr R 
for an expanding, interacting BEC. The presence of this effective mobility edge gives 
rise to two qualitatively different regimes, which might be observed in experiments, 
namely exponential localization for £ in > cr R and algebraic decay of the BEC density 
profile for £ in < a R . 
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Equation (|3"6")) shows that for a speckle potential, the stronger localization is 
obtained for £ in = 3cr R /2 and for this value, L loc increases with a n . It is thus more 
favourable to work with the shortest correlation length of the disordered potential. To 
date, correlation lengths about tr R ~ 0.3/mi have been produced experimentally [48 . 




The first experiments on the expansion of a BEC in a speckle potential [45} |46l [48] 
were operated at £ in ~ cr R /10. In order to reach the regime of Anderson localization, it 
is crucial to lower the interaction energy of the initial condensate in order to increase 
the healing length £ in . One can either use a Feshbach resonance to directly control the 
atom-atom interaction strength or lower the density by lowering the number of atoms 
and/or the radial confinement of the magnetic guide used in Refs. [45, 46, 48 . In the 
latter case, the dynamics of the expanding BEC in the radial direction can play a role. 
As the radial confinement is kept during the expansion of the BEC, we expect the 
radial dynamics to be slow compared to the longitudinal expansion. As a result, the 
radial profile of the BEC would follow adiabatically the local ID (longitudinal) density, 
reducing the dynamics to a quasi-lD problem. However, a part of the potential and 
interaction energies associated to the smooth radial confinement will be converted 
into longitudinal kinetic energy during the expansion. If the BEC is initially in the 
3D Thomas- Fermi regime (/i ^> huj±), the effective de Broglie wavelength of the ID 
expansion can be expected to be A dB — 0.85£ in to be compared to £ in in the pure ID 
case studied here. Hence, for elongated BECs, we do not expect strong differences 
compared to the ID situation we have studied in the present work. A more detailed 
discussion would require further investigation which is beyond the scope of this paper. 

We expect that this work will pave the way for the experimental observation 
of those non-trivial ID localization properties. This work could also be extended to 
higher dimensions (namely 2D and 3D geometries) [98, 99J where similar scenarios can 
be expected but where localization properties would be significantly different. 
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Appendix A. Random and periodic potentials 

In this appendix, we give a couple of details about the inhomogeneous potentials V(z) 
to which the atoms are subjected (see section liOj) i.e. a disordered speckle potential, 
a Gaussian impurity model of disorder and a periodic potential. We plot in Fig. IA1I 
typical realizations of these potentials. 

Generally, we write the potential V(z) as 




V{z) = V R v(z/a R ) 

and the spatial auto-correlation function C(z) = (V(z' + z)V(z')) — (V) 2 as 
C(z) = V£c(z/a K ) 



(A.1) 



(A.2) 




Figure Al. (color online) Typical realizations of the potentials V{z) considered 
in this work, a) Speckle and periodic potentials: the solid line (red online) shows 
a typical realization of a disordered speckle potential with <tr ~ 1.39 X 10 _2 L T f 
and the dotted line (blue online) shows the periodic potential with A = 10 _1 Ltf- 
b) Gaussian impurity model of disorder with w = 2 X 10 — 2 I/tf- 



where V R is the typical amplitude of potential, <t r is the typical space scale (correlation 
length in the case of disordered potentials) and v(u) is a given function characteristic of 
the model of inhomogeneous potential. In this work, we assume that (v(z)) = where 
(.) represents averaging over realizations of disordered potentials or spatial averaging. 

Appendix A.l. Speckle disordered potential 

The main model of disorder we consider is the speckle potential, as it is the 
one used in many experimental studies of disordered Bose-Einstein condensates 
[44l l45l l46l |48l l79l [80] . In brief a speckle pattern is formed by diffraction of a laser 
beam through a rough plate. The intensity of the speckle pattern is proportional to 
the intensity of the incident laser and the correlation function is determined by the 
transmission of the diffusive plate (59] [60] (see Ref. [48] for practical realizations in 
the context of disordered BECs). A disordered potential {e.g. a speckle potential) is 
characterized by its statistical properties, mainly the single-point intensity distribution 
and the two-point correlation function. 

Within the scaling defined above [see Eqs. (|A.ip . (|A.2j) ]. a speckle potential is 
represented by a random function v(u) whose single-point statistical distribution is a 
decaying exponential 

T[v{u)} = exp[-{v(u) + 1}] for v(u) > -1 (A.3) 
and 'P[i>(u)] = otherwise. 

The laser intensity pattern creates an inhomogeneous light shift for the atoms (see 
Fig. lAla ). For a laser which is blue-detuned compared to the atomic resonance line, 
we have V R > and for red detuning V R < [100] . 

For the numerical calculations presented in the paper, we numerically generate a 
ID speckle pattern using a method similar to the one described in Refs. [49] 1101] in 
ID and corresponding to the following reduced correlation function: 

c(u) = sinc(u) 2 . (A. 4) 

Another useful characteristics of the speckle potential in the case of blue detuning 
(Vr > 0) is the average number of peaks with an intensity larger than a given 
value V within a given region of length L TF (see section 1X2]) . Elaborate methods to 
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Figure A2. (color online) Average number of peaks within a range L T f with 
amplitude larger than V for different values of the correlation length of the 
disordered potential and comparison to Eq. (IA.51 ) . 



compute a number of characteristics of speckle potentials can be found in Refs. [59"ll60| . 
Here, we use a simple approximation which suits our purpose. From the probability 
distribution (|A.3|) . we easily find that the probability density that the local disordered 
potential is larger than a given value V is P(V) = exp[—(V+V R )/V R ]. Now, the density 
of peaks (local maxima of the disordered potential) is 1/d where d oc er R is the typical 
distance between two peaks. Therefore, typically, the number of peaks N peaks whithin 

a region of length L TF with intensity larger than V scales as exp V y^' J ■ 
However, this rough estimate does not take into account the interplay between the 
local intensity distribution and the finite correlation length of the disordered potential. 
In the simulated speckle potentials, we find that the typical number of peaks with 
intensity larger than V in a L TF -long region can be approximated by 



N pa ^(V) ~ a I — I exp 



(A.5) 



V R 

with a ~ 0.30 and /3 ~ 0.75, with very good accuracy as shown in Fig. 
Appendix A. 2. Impurity model of disorder 

A model of disorder which is very popular in theoretical studies of quantum disordered 
system is the impurity model [94] : 

3 

where g is a real-valued function peaked at z = 0, of width w and such that 
< g(z) < 1. The locations of the impurities Zj are random and their average 
distance is denoted d. The disordered potential is then formed of a series of impurities, 
all identical but randomly displaced along the z axis (see Fig. lAlb ). Here, we consider 
Gaussian-shaped impurities: 

g(z) - exp(-z 2 /2^ 2 ). (A.7) 
This potential can be realized using ultracold atoms (of another species than the 
BEC) trapped in the Wannier states of the fundamental Bloch band of an optical 
lattice [561 1571 IT02j. 
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From Eq. (| A.6|) , we find that the statistical average of the potential (|A.6[) is 

{V) = T / dx g{x) = ^ Vo © ' (A - 8) 

and the correlation function is 

C ( z ) = % I dx 9(x)g(x + z) = V£c(z/a R ) (A.9) 



with c(u) = y/ir cxp(— x 2 ), V R = \/WVq and cr R = 2w. 

Appendix A. 3. Periodic potential 

For periodic potentials (see Fig. lAla ). we use 

V(z) = V2V R cos(2ttz/A) (A.10) 

corresponding to the mean value (V) = and the standard deviation AV = \V R \. 
Such potentials are currently realized using interference patterns of laser beams in 
geometries with tunable lattice spacings [38| 1103] . 

Appendix B. Phase formalism for the calculation of Lyapunov exponents 
in ID disordered potentials with finite-range correlations 

In this section, we outline the phase formalism method used to calculate the Lyapunov 
exponent (inverse localization length) of a non-interacting particle of energy E = 
ti 2 k 2 /2m in a ID disordered potential with finite-range correlations [M]. The idea 
consists in calculating the propagation of the particle using a perturbation on the phase 
of the wavefunction tp(z). Notice that the perturbation is not on the wavcfunction 
itself. 

The starting point is the Schrodinger equation 

E^(z) = ^^(z) + V(z)^j(z). (B.l) 

Without any loss of generality, we can write the wavefunction and its spatial derivative 
in the form 

${z) = r[z) sin[0(z)] (B.2) 
tp'(z) = kr{z)cos[9{z)\ (B.3) 

where r(z) and 9(z) represent the amplitude and the phase of ip(z), respectively. 
Substituting Eqs. (|B~2|) . (jB~3|) into Eq. (jB~T|) . we find the coupled equations 



e'(z)^k- 2 -^ls^[9(z)} (B.4) 



ln[r(*)/r(0)] = £ dz' sin[20(z')] (B.5) 

Notice that the amplitude r(z) is not involved in Eq. (|B.4[) . It follows that for 
weak disorder [see condition (|B.9|) ]. it can be solved easily in the lowest order of a 
perturbation series of the phase 6(z). We write 0(z) = 9q + kz + 89(z) and we find 

9{z) ~9 a + kz- [ dz' 2m J[ z '^ sin 2 [0 Q + kz n (b.6) 
J h l k 
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in the Born approximation (lower order). Finally, substituting Eq. (|B.6[) into Eq. (|B.5[) . 
we find, in the limit |^| — ^ oo 

ln[r(z)/r(0)] = +l(k)\z\ (B.7) 

where 7(A) = f^^dz C(z) cos(2fcz), where C(z) is the correlation function of 
the disordered potential [see Eq. (|A.2[) ] or equivalcntly 

7 (A) ~ ^ (^j (ka K ) 2 c(2ka R ), (B.8) 

is the Lyapunov exponent. 

Notice that the solution (|B.7[) corresponds to an exponential increase of the 
envelope of the wavefunction ip( z )- I n general, the solution of the Schrodinger 
equation in a disordered potential is the sum of an exponentially increasing function 
and of an exponentially decaying function. In a finite system, boundary conditions 
fix the coefficients and one finds true localized states (i.e. wavefunctions with an 
exponentially decaying envelope for both z — > +00 and z — > —00). In a time-dependent 
propagation scheme, the conservation of probability also imposes that the coefficient of 
the exponentially increasing function vanishes. In the present calculation, since we do 
not impose boundary conditions, only the exponentially increasing function remains 
at infinite distance. 

In spite of this unimportant limitation, this technique provides a very useful 
analytic formula for the Lyapunov exponent, valid for any weak ID disordered 
potential, possibly with finite-range correlations. An important point is that the phase 
formalism technique [M] clarifies the Anderson localization effect in ID disordered 
potentials. Hence, on small distances (say of the order of a R ), the interaction of the 
wavefunction with the disordered potential induces a small perturbation of the phase 
6(z) [see Eq. (|B.4[) ]. but hardly affect the amplitude r(z). Nevertheless, the coupling 
between the phase and the amplitude [see Eq. (|B.5[) ] is crucial and induces at large 
distances (say of the order of L ioc — 1/7) an exponential envelope, characteristic 
of Anderson localization. Finally, the non-interacting localized state of energy E is 
essentially a plane wave of wavenumber k = y^2mE / h 2 modulated by an exponential 
envelope. This makes clear the condition of application of the phase formalism 
approach which requires 7(A) <C A in order for the phase to be only weakly perturbed. 
It follows for Eq. (|B.8|) that this condition reduces to 

V R a R « — (fca R ) 1/2 (B.9) 

TO 

where V n is the amplitude and a R is the correlation length of the disordered potential. 
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